% ========================================================================================================================
% Schniepp Lab, 2018-2021
% Generate a sub peak shape with a mixture of both Lorentzian and Gaussian
% peak features
% ========================================================================================================================

function f=GaussianLorenzNew(AG,AL,gamma,w0,L,w)
%     A. L. Stancik, E. B. Brauns Vibrational Spectroscopy 47(2008)66-69

    LorentzianD=LorentzianDistribution(AL,gamma,w0,w);
    GaussianD=GaussianDistribution(AG,gamma,w0,w);
    f=L*LorentzianD+(1-L)*GaussianD;
end

function f=LorentzianDistribution(A,gamma,w0,w)
    f=(2*A/(pi*gamma))./(1+4*((w-w0)/gamma).^2);
end

function f=GaussianDistribution(A,gamma,w0,w)
    f=(A/gamma)*sqrt(4*log(2)/pi)*exp(-4*log(2)*((w-w0)/gamma).^2);
end